Heavy quarkonia spectra using wave function with gluonic components 
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We calculate the spectra of charmonium and bottomium in an approximation scheme which treats 
hard gluons perturbatively while soft gluons are expanded in a set of localized wave functions. Quark- 
antiquark and quark-antiquark-gluon sectors are included. Reasonable agreement with 2 parameters 
only is found but the spectra are too coulombic. Despite large coupling constant the admixture of 
the quark-antiquark-gluon sector is found to be remarkably small. 
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Quantum Chromodynamics (QCD) is a well established theory describing hadron interactions. Due to asymptotic 
freedom its predictions could be verified in high energy collisions. The bound state problem, however, remains nearly 
intractable. The only ab initio actually available calculations are obtained on the lattice^. There is, of course, a large 
^ body of phenomenological models^ which are quite successful in describing the spectra of hadrons. They are invariably 
based on the assumption that a meson is a quark-antiquark state and a baryon a 3 quark state. Here we present 
an application of a new method applicable to the bound state problem in QCD^. It consists in splitting the gluon 
, momenta into 2 parts: high momenta, treated perturbatively and low momenta, those gluons are expanded not into 
' ■ plane waves but in an orthogonal set of localized functions. As a consequence only a finite number of localized gluons 
,— interact with a heavy quark-antiquark system. High momenta gluons are treated perturbatively within one gluon 
<^ ' exchange approximation. Here we apply this method to the heavy quarkonia (charmonium and bottomium), which 
Oh' allows us to use the non relativistic approximation. To find the eigenstates we postulate a variational wave function 
^ . with 2 components: quark-antiquark and quark-antiquark-gluon. After integrating out gluonic degrees of freedom we 
^ ' obtain a set of two coupled equations describing quark-antiquark system without and in the presence of one gluon. 
, The solution is compared with experimental data. The agreement is not particularly good, our effective potentials 
being too coulombic. However, we use 2 parameters only (quark mass and QCD coupling constant). Moreover the 
I inclusion of quark-antiquark-gluon sector improves the agreement. 
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CSl , II. ZEROTH ORDER HAMILTONIAN IN TEMPORAL GAUGE 
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^""1 . The Lagrange density of QCD is 
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^: 1 - 

g . C = --F^^Fr + ^{iD - m)^. (1) 

^ For the bound state problem in hamiltonian formalism it is convenient to choose a gauge which avoids the problem 
of constraints, the temporal gauge 

^ : ^0-0. (2) 

- - I _J 

In this gauge we have a 3- vector field (A°), its canonically conjugate momenta are simply the components of the 
chromoelectric field 

Af = Ef, (3) 

where the Latin index i takes values 1,2,3 and SU(3) gauge index a takes values 1..8. Introducing chromomagnetic 
field 

Bf = e,,k{d,Al + (4) 
the gluonic part of the Lagrange density is 

^-^-^ (5) 
and the corresponding Hamilton density is simply 

n = ^ + ^. (6) 
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Fields Af and E"" are canonically conjugate and fulfill the commutation relations 

[Anx),E^{y)U=y, = ^S,,S'^%{x~y). (7) 

It should be noted that (6) contains 3 degrees of freedom at each point. The reduction of supernumerary degrees of 
freedom is accomplished by the imposition of Gauss law on the state vectoiii 

G"(x) I * ) = 0. (8) 

For the quarks we use the non relativistic approximation 

H,, = Tr[-^ {pu - 9 A? A? f + J- { p,, + g At )% (9) 
2 TO 2 m 

where A° 2 are the usual Gell-Mann matrices operating on the color index of the first (second) particle. 



III. HIGH-LOW DECOMPOSITION 



In order to reduce the problem to a numerically manageable form we have to split our Hamiltonian eqs (6) (9) into 
a zeroth order solvable Hamiltonian and the rest to be treated perturbatively. We split the gluons into 2 groups: high 
momentum (hard) gluons (with pi > A), where A is the cut off parameter, and low momentum (soft) gluons. In high 
momentum sector we use the conventional approach - the interactions are treated as a perturbation. In the lowest 
order we have free gluons and quarks interacting via truncated coulomb potential^. Notice that we have dropped 
spin-spin and spin-orbit terms. 

Vc = ,^ ""^,2 ' I Pi - P2|>A, (10) 

where pi, p2 are the momenta of quark and antiquark respectively and as = |a. For soft gluons we use a different 
approach. Instead of expanding the field into plane waves, we expand them into a complete set of localized functions 
in position space. We have chosen to use the Fourier expansion in a cube —A > > A in momentum space. In 
position space our wave functions are the product of 3 well known sinc{x) = sin(x)/x functions 

es(a;, y, z) = ( — ) 2 sincinxTT — xK) sinc{nyT: — yA) sinc{nzTr — zA), (11) 

where ft = (nx,ny,n2),nx,ny,nz are integers (positive and negative). These functions are distributed on a lattice 
with lattice spacing tt/A. They overlap very little and, if A is small enough, the width of the peak will be much larger 
than the rms radius of heavy quark-antiquark system. We have found that A= 0.15 - 0.3 GeV was sufficient and the 
results vary very little with A. Thus for soft gluons we have the expansion 

Anr)^yenmi^ (12) 



and similarly for the chromoelectric field 



Enr)^yen{r)Elf, (13) 



where the Af and £'f are the operators of localized gluons and fulfill 

[AU , 4™] = iS,,S'^''6:,{n-m) (14) 

Since the wave function of quark-antiquark system " fits" within the wave function at one lattice site (chosen to have 
fi — (0,0,0) all others gluons will interact very little with quark-antiquark system. Thus it is reasonable to take into 
account, in the zeroth approximation, only the gluons with n = (0,0,0). In fact the quark-antiquark wave function 
is so narrow that we can approximate e,-j(r) by its value at the origin. Our zeroth order hamiltonian describes a 
quark-antiquark pair with a truncated Coulomb potential interacting with gluons localized at the origin 

Ho^—{pi + g{-)--A,f+—{p2-a{-)-^A2f + V,+^-^ + 15 
2m7r 2TO7r 2 2 
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where we have used 

XIA^ = ii (16) 

and similarly for the antiquark. We have also dropped the index n = (0,0,0). Bf is the chomomagnetic field with 
derivative terms dropped (functions eij(r) are slowly varying within the range of quark-antiquark wave function) 

Bt = g Sijk r"' A". Al (17) 

jnhc .jj.^ structure constants of the SU(3) group. We have reduced the problem to a quark-antiquark system 
interacting with 3 x 8=24 gluons localized at the origin. This is numerically tractable by a method widely used in 
atomic and nuclear physics: Ritz variational method. Not all eigenstates of this hamiltonian are physically acceptable 
- the wave function must fulfill the Gauss law eq. (8), i.e. the state must be colorless. After elimination of center of 
mass motion we get finally 

where p is the relative momentum of the quark-antiquark. 

IV. VARIATIONAL ANSATZ 

The mass of the quark-antiquark system is the energy difference between the energy eigenvalue of Hamiltonian 
Ho and the energy of vacuum fluctuations i.e. the lowest energy of gluonic part of Hq. We have chosen to work in 
Schrodinger representation with 

Et = -^doAt, (19) 
For the vacuum part of the problem we have chosen a variational gaussian Ansatz 

^yac{A) = Nexp{ -pyac ^2°) 

(3 being the variational parameter. We obtain 

4 2^a^A 

Pvac = , = 0.168 GeV, (21) 

(37r) 3 

12 63 ah A 

= -^-^ = 1.511 GeV. (22) 

773 

For the quark-antiquark system we postulate a 2 component wave function 

*(p,i) = 6ij Mp) exp{ -/? ^^i^} +f2{p)^ exp{ - /3 Z^Ifilil} (23) 

The first part describes a quark antiquark system with no gluons (gaussian i.e. lowest energy state of the oscillator) 
while the second part represents a quark-antiquark-gluon system (the gluonic part is proportional to the first excited 
state of an oscillator). Both parts arc colorless. The variation is to be done with respect to /3, fi{p) and /2(p) with 

the condition that the norm of the wave function \l/( p , A ) be one. After solving the variational problem for fi{p) 
and /2(p) and integrating out gluonic degrees of freedom we have a system of two coupled equations 



where 
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~ 4 32 

hii - ^ + V.{r) + %f- + 6 /? + '-(^ (25) 

31 5^ A3 13/3 112^2 A3 

/.I. = (-^) P (27) 
3 p /I 7r2 

A - 

h2i = (^)P- (28) 

Notice that diagonal terms differ by a constant only and that the coupling terms are proportional to p. This set of 
equations was solved numerically by expanding fi{p) and f2ip) in the basis of Laguerre polynomials. We found that 
a basis of 16 polynomials was sufficient. 



V. RESULTS AND DISCUSSION 



We have calculated spectra of charmonium and bottomium. In both cases we used 

A = 0.15 GeV (29) 

In our ab initio calculation we have 2 free parameters only: the mass of the quark and the coupling constant. For 
charmonium we found 

rric = 1.8604 GeV, (30) 
= 0.973, (31) 

by fitting the 2 lowest L=0 states. Since we have no spin-spin interaction in our non relativistic approximation ^So 
and ^Si are degenerate, we used weighted the average^. All the experimental data are from refi^-. 
The results are summarized in Table 1 and Table 2 and shown in Fig. 1 and Fig. 2. We show the results for a 
calculation with quark-antiquark sector only {Sqq), with quark-antiquark-gluon only (Sqqg) and full coupled channels 
(Sfuu)- We consider the agreement reasonable considering that only 2 parameters were used. In our approximation the 
system is too coulombic but the inclusion of quark-antiquark-gluon sector improves the agreement with experimental 
data. To our surprise the admixture of quark-antiquark-gluon sector is very small; for the ground state we have 

Pqqg = 0.0017 (32) 

despite using a large value for a^c ~ 1- This could be an explanation why phenomenological models with constituent 
quarks only are so successful^t. For bottomium the results are similar. We find 

nib = 5.1896 GeV, (33) 
a^b = 0.66114. (34) 

Notice that the coupling constant is 30% lower. Again the spectrum is too close to a purely coulombic one but 
the inclusion of quark-antiquark-gluon sector improves the agreement (Table 3-4 and Fig 3-4). The admixture of 
quark-antiquark-gluon sector is very small 

Pqqg = 0.0004. (35) 

The main feature of this calculation is the surpisingly small admixture of quark-antiquark-gluon sector. However 
the long range part of the effective potential is clearly wrong. The obvious improvement would be to take into account 
space variation of gluon wave function and the inclusion of gluons located on neighbour sites on the lattice. While 
the computation of spectra would be vastly more complicated it is certainly feasible. In conclusion it is seen that it 
is possible to perform bound state calculations within QCD and improve them step by step. 
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Table 1 : S levels in cc (in GeV) 
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Table 3 : S levels in bb (in GeV) 
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Table 2 : P levels in cc (in GeV) 
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Table 4 : P levels in bb (in GcV) 
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